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Abstract 



We demonstrate the consistency at the next-to-leading-logarithmic (NLL) level of a factorization theorem based on Soft- 
Collinear Effective Theory (SCET) for jet shapes in e + e~ collisions. We consider measuring jet observables in exclusive multijet 
final states defined with cone and kj-type jet algorithms. Consistency of the factorization theorem requires that the renormalization 
group evolution of hard, jet, and soft functions is such that the physical cross section is independent of the factorization scale fi. 
The anomalous dimensions of the various factorized pieces, however, depend on the color representation of jets, choice of jet ob- 
servable, the number of jets whose shapes are measured, and the jet algorithm, making it highly nontrivial to satisfy the consistency 
condition. We demonstrate the intricate cancellations between anomalous dimensions that occur at the NLL level, so that, up to 
power corrections that we identify, our factorization of the jet shape distributions is consistent for any number of quark and gluon 
^ jets, for any number of jets whose shapes are measured or unmeasured, for any angular size R of the jets, and for any of the algo- 
— j rithms we consider. Corrections to these results are suppressed by the SCET expansion parameter A (the ratio of soft to collinear or 
collinear to hard scales) and in the jet separation measure 1/f 2 = tan 2 (7?/2)/ tan 2 (i^/2), where iff is the angular separation between 
jets. Our results can be used to calculate a wide variety of jet observables in multijet final states to NLL accuracy. 
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Final states that contain several jets are important Standard 
Model backgrounds to many new physics processes in high- 
energy colliders, in addition to serving as sensitive probes of 
Quantum Chromodynamics (QCD) itself over a wide range of 
energy scales. The structure of jet-like final states contains sig- 
natures of the hard scattering of parton-like degrees of free- 
dom, the branching and showering at ever lower energies, and 
hadronization at the lowest scale Aqcd- Probing the structure 
of jets both teaches us about QCD and can help us to distin- 
guish jets of Standard Model origin from those that are truly 
signatures for new physics. 

The presence of multiple scales governing jets is at once 
an opportunity to probe many aspects of their physics and also 
a challenge due to the generation of large logarithms of ratios 
of these scales spoiling the behavior of perturbation theory. A 
powerful framework to separate physics at different scales and 
to improve the behavior of perturbation series is effective field 
theory (EFT). EFTs aid in factorizing an observable dependent 
on multiple scales into pieces each sensitive to a single energy 
scale. Renormalization group (RG) evolution of these pieces in 
EFT achieves resummation of large logarithms to all orders in 
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perturbation theory. Factorization also allows the disentangling 
of perturbative and non-perturbative physics |Q]|2). 

Soft-Collinear Effective Theory (SCET) 11 15) has had 
considerable success in applications to many hard-scattering 
cross sections |7| and jet cross sections. SCET separates de- 
grees of freedom in QCD into distinct soft and collinear modes, 
expanding the full theory in a parameter A that characterizes 
the size of collinear momenta transverse to the jet direction, 
and provides a framework to factorize cross sections into sepa- 
rate pieces coming from interactions at hard, collinear, and soft 
scales. This was done in SCET for event shape variables us- 
ing hemisphere jet algorithms in <? + e~ colliders |8] and for 
"isolated Drell-Yan" (where central jets are vetoed) in hadron 
colliders [ 10 1 . In addition, there has been progress in under- 
standing how to implement jet algorithms other than the sim- 
ple hemisphere jet algorithm in SCET. In ATI [12 1. total two-jet 
rates where the jets are defined by Sterman-Weinberg jet algo- 
rithms were computed at NLO. These results were extended to 
the cases of the exclusive k T and JADE algorithms in fOl . 

In most applications of SCET to exclusive jet cross sections 
considered to date, there are two back-to-back jets. (Recently 
Ref. [14 1 considered direct photon production in hadron col- 
lisions, involving three collinear directions.) In this work we 
consider for the first time exclusive N-jet final states with arbi- 
trary N >2 for the SISCone 1 15 1, Snowmass 1 16 1, inclusive k T 
1 17 1, anti-kr fHfl . and Cambridge- Aachen |19| jet algorithms. 
We find that a new feature that arises when more than two jets 
are present is that the parameter A is not in itself sufficient to 
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ensure factorization. In particular, factorization is valid to lead- 
ing order in A and in a jet separation measure l/t, where t is 
defined by 

t = tanQfr/2) 
tan(/?/2) ' 

with R the angular size of a jet as defined by a jet algorithm 
and iff the minimum angle between two jets. This is due to 
the fact that jets need to be both well-collimated (i « 1) and 
well-separated (f » 1). The latter requirement is trivial for 
back-to-back jets since l/t = for if/ = n. 

Our analysis applies not only to the total A^-jet cross section, 
but also in the case that jet observables are measured on some 
number M < N of the jets. We will illustrate the measurement 
of angularities T a (cf. j20l|2T|), defined by 



J ieJ 



(2) 



where Ej is the energy of the jet J, the sum is over particles 
i in the jet, and p' T and 77, are the transverse momentum and 
(pseudo-)rapidity of particle i with respect to the jet axis. How- 
ever, most of our results do not depend on this choice of observ- 
able, and we organize the calculation such that other observ- 
ables can be easily implemented. Distributions of jet shapes 
such as angularities contain logarithms of r fl that become large 
in the limit r fl — > 0. The factorization theorem we present pro- 
vides the basis for resummation of these logarithms to all orders 
in perturbation theory. 

Factorization of event shape distributions in SCET was prov- 
en in Refs. Il22ll23l . and factorization for multijet observables 
defined with arbitrary algorithms was considered in Ref. [24|. 
The extension to the more general case that we consider in- 
volves the straightforward combination of the techniques devel- 
oped in these papers and will be derived in detail in Ref. 1 25 1 . 
In this work we demonstrate that, after intricate cancellations 
among the various contributions to the jet and soft functions, 
consistency of the factorization theorem is satisfied at next-to- 
leading logarithmic (NLL) accuracy. In order for the factoriza- 
tion theorem to be consistent, the hard, jet, and soft functions 
defined must satisfy a strong condition on their anomalous di- 
mensions: 
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for any number N of total jets and M of measured jets, and any 
color representation of each jet. This consistency condition is 
made even more nontrivial by the potential dependence of the 
jet and soft anomalous dimensions on the jet algorithm parame- 
ters. In this Letter we demonstrate that Eq. ([3]) does in fact hold 
for arbitrary numbers, types, and sizes of jets in the final state, 
up to certain power corrections we are able to identify. 

We begin in Sec. |2]by defining the phase space cuts needed 
to implement our choice of jet algorithms. In Sec. [3] we then 



present the factorization theorem for A^-jet events and define the 
hard, jet, and soft functions, and identify power corrections to 
the factorization. In Sec.|4]we give the form of the RG evolution 
equations obeyed by the factorized functions. In Sec.[5]we sum- 
marize the results of all the anomalous dimensions needed for 
NLL running and demonstrate how they intricately satisfy the 
consistency condition Eq. This requires calculating only 
the infinite parts of the bare functions. We give the finite pieces 
of the jet and soft functions (which are not needed at NLL) in 
Ref. [25]. In Sec.[6]as an example we calculate quark and gluon 
angularity jet shapes in 3-jet final states with logarithms of r a 
resummed to NLL accuracy. 

2. Phase Space Cuts and the Jet Algorithm 

Two general categories of jet algorithms, cone algorithms 
and recombination (kj-type) algorithms, are commonly used to 
find jets. For a jet composed of two particles, as in a next-to- 
leading order description, the phase space constraints implied 
by each type of algorithm become very simple. In this work 
we deal with the common forms of cone and (inclusive) kx- 
type algorithms; our cone algorithms include the Snowmass and 
SISCone algorithms, and our recombination algorithms include 
the inclusive k T , Cambridge-Aachen, and anti-k T algorithms. 
Cone algorithms require each particle to be within an angle R 
of the jet axis, while recombination algorithms require the angle 
between the two particles to be within an angle D of each other. 
If we label the jet axis as n and its constituent particles as 1 and 
2, then the algorithm constraints for a two-particle jet are: 



cone type: #i„ < R and 6o n < R, 
k T type: Qyi < D . 



(4) 



For the parts of the jet and soft functions that we give in this 
work, we find that the functional form is the same for cone-type 
and kj-type algorithms in terms of the angular parameter R or 
D. Therefore, we will use the more common R in writing down 
the jet and soft functions, but we note here that the functional 
form is the same for kr with the replacement R — > D. 

Note that, while all algorithms that we consider fall into 
one of the two constraints in Eq. Q at NLO, at higher orders 
the various algorithms will behave differently. Without taking 
this into account, we have no guarantee that we can resum all 
logarithms of jet algorithm parameters correctly. This is not a 
problem we solve in this paper. In this paper, we resum loga- 
rithms of jet observables in the presence of phase space cuts due 
to an algorithm, demonstrate that the factorization theorem and 
NLL running are valid and consistent, and identify the power 
corrections to this statement. 

At the hard scale, we match an N-leg amplitude in QCD 
onto an A^-jet operator in SCET, meaning we must enforce that 
the number of jets is fixed to be N. To enforce that we have 
no more than N jets, we require that the total energy of parti- 
cles that do not enter jets to be less than a cutoff A. To enforce 
that we have at least N jets, we need that pairwise each jet is 
well separated from every other jet. The requirement of consis- 
tency of NLL running will give a quantitative measure of this 
separation requiring that t » 1 . 



3. Factorized Jet Shapes in A 7 - Jet Production 

The cross section for e + e~ annihilation to N jets at center- 
of-mass energy Q, differential in the jet three-momenta P,- of 
the jets and in the shapes of M of these jets, is given in QCD by 

dcr 



(5) 



M N 
X &n(3(X))-N f[ S(T' a - T a (Ji)) [J 6\Pj - P(Jj)) , 

i=\ j=\ 

where 7,- is the ith jet in X identified by the jet algorithm J~. The 
Kronecker delta restricts the sum over states to those that are 
identified as having N jets by the algorithm. The final state is 
produced by the QCD current f = qy^q, and is the leptonic 
part of the amplitude for e + <?~ — > y* . 

To factorize the cross section Eq. d3), we begin by matching 
the QCD current f onto a set of A^-jet operators in SCET. These 
operators are built from quark and gluon jet fields, 



X n = W% , &t = - Wl(P ± + A^)W n , 



(6) 



where £„, A n are collinear quark and gluon fields in SCET, and 
W„ is a Wilson line of the 0(1) component h ■ A„ of collinear 
gluons, 



W„(x) = ^ exp 



perms 



n ■ A n (x) 



(7) 



We have made use of the label operator 7* which picks out the 
large 0(1) n-p and 0(A) p ± components of the label momentum 
p of collinear field in SCET. We will not need to construct the 
A^-jet operators explicitly, but bases of 2, 3,4 jet operators have 
been given in [7, 26, 27], respectively. 

To describe an A^-jet cross section, we construct an effec- 
tive theory Lagrangian by adding N copies of the collinear La- 
grangian in SCET (in N different light-cone directions n,-) to- 
gether with one soft Lagrangian. In each collinear sector, we 
redefine collinear fields by multiplying by Wilson lines of soft 
gluons to eliminate the coupling of soft gluons to collinear modes 



in the leading-order SCET Lagrangian [6|, g n = Y„^ n 
A n = J/„A,f , where 



(0) 



Y n (x) = P exp 
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dsn ■ A s (ns + x) 



and 



(8) 



with A s in the fundamental representation, and J/ similarly de- 
fined but in the adjoint representation. 

Performing the above steps in Eq. |5]) for the jet shape dis- 
tribution, the details of which we report in [25 1, we obtain the 
factorized form in SCET, 



dcr 



daW 



-ff(p 1 ,...,p JV )f]4^ 



M „ 



where cr (l)) is the Born cross section for e + e~ — > A^ partons, 
H = 1 + 0(a s ) is the hard coefficient given by the matching 
coefficient of the SCET A^-jet operator, and J and S are jet and 
soft functions. The superscripts f; denote the color representa- 
tion (corresponding to a quark, antiquark, or gluon) of the jet 
corresponding to the /th leg in the A^-jet operator. We number 
the legs so that i = 1, . . . , M are the jets whose shapes we mea- 
sure, and the remainder j = M + 1, . . . , N are left unmeasured. 

The quark and gluon jet functions for jets whose shapes are 
measured are defined by 

■C(t/) = ^Tt2 f d\e-^h nmXn)yl 

Nc l^J 2n J 2 (10a) 

X (0\XnAx) \*n) (X n \Xn,ca(0) |0> 6(tj - T a (J(X n ))) , 
<jj 



(J(X„))-l 



x — <0| gBZ(x) \X„) (X n \gB^(0) |0) 6(tj - T a (J(X n ))) , 

(10b) 

where the traces are over color and spinor indices, and D is 
the number of dimensions. The sums are over states in the n- 
collinear sector. The label direction and energy n, a> are cho- 
sen to match the jet momentum P. We have factored the Kro- 
necker delta in the full cross section Eq. |5]) restricting the sum 
over states to those with Af jets according to the algorithm J' 
into individual restrictions that there is precisely one jet in each 
collinear sector. The delta functions of tj restrict the angularity 

of the jet J identified in the state X n by the jet algorithm. The 

f 

jet functions /„ \ b) . for jets whose shapes are left unmeasured 
are given by Eq. (|T0j> without the delta functions of tj. 

The soft function, meanwhile, is given by matrix elements 
of Af soft Wilson lines in each of the collinear directions n, and 
color representations r, of the ith jet. For arbitrary Af, multiple 
color structures may appear, and if so there is an implicit sum 
over multiple hard functions H and soft functions S in Eq. 
An A^-jet soft function takes the general form, 



1 

Sn^)) = - J] 6 nmXs)) Y] 5^ - T l a (X s )) 



(11) 



\YZ r ---Y^(0)\X s )(X s \Y r n \---Y^(0)\0) 



(9) 



where N normalizes the soft function to 6(t^) ■ ■ ■ 5(t^) at tree 
level. There is an implicit contraction of color indices which 
we have left unspecified. The whole soft function is color sin- 
glet. Note that the sum over soft states is restricted so that soft 
particles do not create an additional jet when the jet algorithm 
is run on X s . r' a (X s ) is the contribution to the jet shape from soft 
particles which are actually in the jet 

The factorization of the cross section Eq. |9| is valid in the 
following limits of QCD: 

1 . The SCET expansion parameter A, determined either by 
the jet shape r a for measured jets or the jet radius R for 
unmeasured jets, must be small. In other words, each jet 
must be well collimated. 



2. The separation between any pair of jets must be large. We 
will find that the natural measure for this separation is the 
variable t = tan(i^/2)/ tan(R/2), where ijj is the minimum 
angle between two jet directions, t must be large, that 
is, jets must be well separated in order for us to factor 
the N-jet condition in the full cross section Eq. Q into 
N individual 1-jet conditions in each collinear sector as 
in Eq. ( fT0| ) and a no-jet condition in the soft sector as in 
Eq. ( flT) . This approximation is inevitable because each 
jet function 7, already approximates all radiation emitted 
by other jets as coming from a Wilson line W ni along the 
exactly back-to-back direction whereas the hard and 
soft functions know the directions of all N jets exactly. 

3. The energy of all particles not included in a jet must be 
of the order of soft momenta. This is so that setting the 
label energy on each of the jet fields in Eq. ( fT0) > to be 
equal to the total jet energy is correct at leading order in 
A. In particular, the energy cut parameter A on energy 
outside of all jets is required to be soft, A ~ A 2 Ej. 

4. Power corrections associated with the jet algorithm are 
small. For instance, setting the jet axis equal to the la- 
bel direction n is valid up to 0(A 2 ) corrections, which in- 
duce corrections to the jet shape T J a which are subleading 
for a < 1 1 20 , 2~3l l28l . Similarly, assuming soft parti- 
cles know only about the total collinear jet momentum 
by the time they are included or excluded from a jet in- 
duces power corrections to t ( { that are power suppressed 
for sufficiently large R. 

We go into greater detail about these approximations in (25]. 

4. Renormalization Group Evolution 

The functions that we consider either renormalize multi- 
plicatively or through convolutions in r. The multiplicative 
form of a renormalization group equation (RGE) obeyed by a 
function F is 

fi^-F(M)=y F Qi)F(ji), (12) 
dp 

with the anomalous dimension of the form 



y F (p) = r F [a]ln 



OJ 2 



■y F [a] . 



This RGE has the solution 

FQx) = U F (p,p )FQi ) 

where 

U F ( M , M0 )=e K ^(^ ) j 



(13) 



(14) 



to F Qi,fto) 



(15) 

where we define a> F , K F below in Eq. ( pO) . The convolved form 
of an RGE obeyed by functions F that depend on the observable 
is 

d r 

(16) 



fi d^p F ^ T, ^ > = f dT'y F (T-T';p)F(T';p), 



where 



y, ( r: u) = j r>M IrA + y F [a]) 6(t)- — r F [a] 

U 2 I JF 



0(T) 



(17) 



The solution to this RGE is ll8ll29ll30ll3T1[32l 

F(r,/i) = J d T 'U F (T-T';p,po)F(T';p ), (18) 



where 



p K F +y E ai F . \j F aj F 



(19) 



We note that the anomalous dimensions y F (p) and y F (r; p) in 
general also depend on the jet algorithm parameters R and A 
which we have made implicit. 

The part of the anomalous dimensions in Eqs. ( p"3j ) and fP7) i 
multiplying ln(p 2 /u> 2 ) is proportional, to all orders in a s , to the 
cusp anomalous dimension T(a s ), given to 0(a s ) by r(a s ) = 
ajn. With one-loop results for the anomalous dimensions, and 
using the two-loop form of the cusp anomalous dimension, the 
RGE solutions are accurate to NLL order. In Eqs. ( fT5] l and (19) , 
(Op, K F are given by 



2 r as(M) da 
u F (pi,p ) = — ■ 7 ^T F [a] 



K F (p 



,Mo) = I 

Ja s (b 



■"^ da 
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da 



y F [a] 



+ 2 



r a " w da „ r r° 

Ja s (ji ) Pl a l Ja s (ji 



da 



(20a) 



(20b) 



where f3[a] is the beta function of QCD. We define j F = 1 for 
RGEs of the form Eq. (TJ). 

We will find that the hard function can be written as a sum 
over functions that each obey a multiplicative renormalization 
group equation. The unmeasured jet function also obeys a mul- 
tiplicative RGE, while the measured jet function obeys a RGE 
with a convolution over r. The soft function, whose struc- 
ture we will discuss in detail, can be decomposed into terms 
which obey multiplicative RGEs and terms which obey con- 
volved RGEs. 

In the next section we outline the calculations necessary to 
obtain all the above anomalous dimensions to 0(a s ). 

5. Anomalous Dimensions and Consistency of Factorization 

In this section we discuss the calculation of the one-loop 
hard, jet, and soft anomalous dimensions and the form of the 
anomalous dimensions in Tableland demonstrate that the con- 
sistency condition, Eq. is satisfied to one-loop order, to 
leading order in the approximations we enumerated above. This 
is already an intricate test whose satisfaction turns out to be 
highly nontrivial. Having verified this condition, we proceed at 
the end of the Letter to give an application of NLL resumma- 
tion of the jet shape distribution made possible by our one-loop 
calculation of the anomalous dimensions. 

5.1. Hard Function 

The hard function H in the factorized cross section Eq. (|9| 
is given by the square of the Wilson coefficient in the matching 



of the Af-parton amplitude in QCD onto an Af-jet operator in 
SCET, 

(N\qTq\0) = (N\C N O N \0) , (21) 

where the right-hand side is actually a sum over many possi- 
ble N-jet operators built from the jet fields in Eq. |6) and soft 
Wilson lines Eq. (jH). The allowed basis of operators On is de- 
termined by gauge and Lorentz symmetry. If there is only one 
operator, the hard function is simply H = \Cn\ 2 - 

The one-loop anomalous dimension of the N-jet matching 
coefficient Cn can be determined from calculations existing in 
the literature, for example, Table III of Ref. [33 1. For an opera- 
tor with N legs with color charges T,, the anomalous dimension 
of the matching coefficient Cn is 



Tfrc^in^ + iy/te) 

U>i i 



7c N (a s ) =~Yj 
;=i 

1 v-^ i-n.i-n.i- ;0 4 
-2 r(as) Z T <' T ; ln { 

where y, is given to 0{a s ) for quarks and gluons by 
3a s C F a s UC A -4T R n f 



(22) 



7q 



In 



Yg 



(23) 



The anomalous dimension of the hard function itself is then 
given by yn = jc N + 7c an d can be written as 



(24) 



!=1 



Because the hard function obeys a multiplicative RGE, each 
term in the hard function obeys a multiplicative RGE, and so 
each term in Eq. p4| ) has the form Eq. ( |13) . Each H' has a> = 
(jj h while Y[a] = for IP™, as listed in Table[l] 

5.2. Jet Functions 

The quark and gluon jet functions are given by Eqs. ( |10a| i 
and (JTObJi and are calculated from cutting all possible diagrams 
at a given order in a s correcting a collinear propagator with 
label momentum ton. The jet functions include phase space 
restrictions on the final-state particles from the cut requiring 
that only one jet is produced. When we cut through a single 
propagator, the solitary parton in the final state is automatically 
in the jet, but these diagrams turn out to be scaleless and thus 
zero in dimensional regularization. For the cuts through loops, 
two collinear particles are created in the final state, and both 
particles are in the jet if Eq. Q is satisfied. If Eq. Q is not 
satisfied, we require one of the particles to have energy E < A, 
so that only one jet is produced by the final state. Additionally, 
for jets whose shapes are measured, we include a delta function, 
6(t j — T a (J(X))), measuring the jet shape for the particles in the 
jet. The restrictions on unmeasured jet functions are the same 
as the measured jets except for this delta function. 

We report here the results of calculating only the infinite 
parts of the relevant loop graphs in dimensional regularization, 
in D — 4 - 2e dimensions, in the MS scheme. We give the finite 
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Table 1: Anomalous dimensions of hard, jet, and soft functions. The cusp 
parts Tf and non-cusp parts yp of the anomalous dimensions for hard, unmea- 
sured jet, measured jet, and soft functions are given, along with the constant 
]F appearing in Eqs. (17) and (20a) . T is the cusp anomalous dimension, given 
to one-loop by T = a s jn. The pieces y, for quarks and gluons are given by 
Eq. (23). The three rows for the soft anomalous dimensions are organized to 
correspond to the three groups of evolution factors given in Eq. (32) and are 
given in the limit 1/r 2 — > 0. 



parts in 1251 . Our calculations give anomalous dimensions for 
quark and gluon jets y 1 , of the form Eq. ([13) for unmeasured 
jets and y^ija) of the form Eq. fP7] > for measured jets, with the 
values given in Table [T] 

In the measured jet function, we find that the zero-bin sub- 
traction plays a key role. The zero-bin subtraction removes 
doubly-counted regions of phase space from the "naive" con- 
tributions to the jet function |34|. For the measured jet func- 
tions, the naive contributions to the anomalous dimension only 
depend on 5(j a ) and do not contain (1 /t cj ) + distributions. How- 
ever, the zero-bin contribution to the anomalous dimension con- 
tains non-trivial r a dependence away from r a = 0, and it is only 
by performing the zero-bin subtraction that we obtain the cor- 
rect running of the measured jet function. 

When the final-state particles in the jet function do not pass 
the cuts in Eq. Q, only one particle is in a jet. In this case the 
contribution to the jet function is power suppressed by 0(A/u>), 
since a collinear parton must have E < A to be outside of the 
jet. This power contribution is not power suppressed in the 
naive contribution alone, but only after the zero-bin subtraction. 
Additionally, the zero-bin removes the dependence of the mea- 
sured jet function anomalous dimension on the jet algorithm 
parameter R. For unmeasured jets, the zero-bin is a scaleless 
integral, and the R dependence remains in the unmeasured jet 
function. 

Tabulating the results, we find the anomalous dimensions 

are 



jj. = F(q-. s )T 2 In 



Loj tan 2 f 



(25) 



for unmeasured jet functions, and 



= T 2 



2F(a. s )T 2 



+ Ji 



l-a 



~ 0(T a ) 



(26) 



for measured jet functions. 
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5.3. Soft Function 

The soft function in an A^-jet cross section is given by Eq. (jTTJ, 
containing matrix elements of N soft Wilson lines in the N jet 
directions, with each Wilson line in the color representation of 
the corresponding jet. At 0(a s ), this soft function is given by 
a sum over cut diagrams represented in Fig. [T] The blob rep- 
resents the jet in direction and we leave implicit the phase 
space cuts needed for each diagram. We use Feynman gauge, in 
which each diagram is proportional to n, • rij. (Note this allows 
us to drop graphs with i = j or i = k since n 2 = 0.) 

To calculate the soft function, we must implement phase 
space cuts on the soft gluon in the final state requiring that it 
either be in a jet or not produce a new jet (i.e., it has energy 
less than A). The soft function is a sum over contributions from 
all pairs of directions i and j that exchange the soft gluon, and 
we calculate the total contribution with i and j fixed before sum- 
ming over directions. A natural way to organize the phase space 
of the soft gluon in the final state is as follows: 

(1) The gluon enters a measured jet and contributes to t*(X s ). 

(2) The gluon enters an unmeasured jet and has any energy. 

(3) The gluon is not in any jet and has energy E < A. 

We name contribution (1) 5™ eas (r^), where the subscript ij de- 
notes that the gluon goes from i to j. Regions (2) and (3) do not 
contribute to the angularity of any jet and just give an additive 
contribution s™ n - meas to the coefficient of 5{r\) ■ ■ ■ <5(rf ) in the 
full soft function S(j\, . . . , t%). Contribution (3), however, is 
very awkward to calculate, as we must integrate over a phase 
space with many "holes" (corresponding to the jets) removed, 
resembling Swiss cheese. It is easier to reorganize contributions 
(2) and (3) into the following form: 

(A) S 1 ™ 01 : the gluon is anywhere with energy E < A. 

(B) S 1 *.: the gluon is in jet k with energy E > A. 

(C) Syi the gluon is in jet k with energy E < A. 

Then, the unmeasured soft gluon contribution S ™ meas (the sum 
of (2) and (3) in the original list) is given by the combination 



c incl , 



(27) 



k=M+l 



k=\ 



In the first term, coming from region (A), we filled in the holes 
in the Swiss cheese-like region (3) in the original list, allowing 
the soft gluon to go anywhere with energy E < A. We compen- 
sated by adding the second term given by region (B) containing 
gluons with energy E > A inside unmeasured jets (part of the 
original region (2)) and subtracting the third term from region 




Figure 1: Soft Function Diagrams. A gluon exchanged between jets i and j 
crosses the cut which imposes phase space restrictions due to the jet algorithm. 
The blob represents the jet in direction k, which the gluon may enter or not. 



(C), removing gluons with E < A inside measured jets, which 
are already correctly accounted for in 5™ eas (r*). 
The total soft function at 0(a s ) is then given by 



M 



M 



k=\ l=\ 



l=\ 



(28) 



unmeas 



{\8{T l a ) 



Note that the second line is independent of the jet shape. This 
contribution is universal and will appear in any A^-jet cross sec- 
tion in which some of the jets defined by a particular jet algo- 
rithm are not measured. 

The contributions of the measured jet piece 5™ eas (r*) to the 
anomalous dimension of the soft function are given in Table [2] 
separately in the cases that k = i or j and k + i, j. These contri- 
butions are given by the form Eq. ( fP7j ), with the values given in 
Table [2] The results are given in terms of the distance measure 
tjj = tan(i^,j/2)/ tan(7?/2) between jets of size R separated by 
an angle ^y, and the angle jS, ; between the ik and jk planes. For 
well-separated jets, the contributions to the non-cusp part of the 
anomalous dimension are suppressed by 1/f 2 . 

The "inclusive" contribution S„ for a soft gluon going 
anywhere with energy E < A contributes a term to the soft 
anomalous dimension given by the general form Eq. ( |13) , with 
values given in Table [2] 

Finally, for the contributions of soft gluons entering jets 
with E > A or E < A in (B) and (C) in the list above, we 
can combine the last two terms in Eq. |27} using the following 
observation. The sumSy+Sy is the contribution of a soft gluon 
entering jet k with any energy. The phase space integral for this 
contribution contains a scaleless integral (of energy from to 
oo), and so this sum is zero in pure dimensional regularization. 
Thus we can set S'jj = S*p and the last two terms in Eq. ( |2"7] > 
add up to the contribution of a soft gluon entering any jet with 
energy E > A. These contributions can again be split up into 
those with k — i or j and k + i, j. They contribute parts to the 
soft anomalous dimension falling into the form Eq. ( fT3"j ), with 
values in Table [2] The non-cusp pieces are again suppressed by 
1/f 2 for well-separated jets. 

Using the contributions described above, we sum over di 
rections i and j and obtain the anomalous dimensions for S 
and 5 unmeas ? which we record in Table [2] 

The soft function obeys the renormalization group equation 

ft-^S(Tu...,T M ;M) = J dT\---dT' M S(T\,...,T' M \n) 

X7s(ti -t\,...,t m -t' m ;h). 

Because the soft function at 0(a s ) in Eq. ( |28] i is a sum of terms 
that depend non-trivially on at most one jet shape, the anoma- 
lous dimension can be decomposed as 

ys(Ti, . . . , t m ; n) = y—(M) Sin) ■ ■ ■ 5{t m ) 

M M 

+Eto™oI~K>' (30) 
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Table 2: Soft Anomalous Dimensions. Contributions to the anomalous dimen- 
sion of the soft function are given for soft gluons emitted by jet i or j and 
entering jet k (with k — i or j in the first row and k # £, / in the second) and 
being measured with angularity r^; soft gluons emitted by jet i or j in any di- 
rection with energy E < A in the third row; and soft gluons emitted by jet i or j 
and entering jet k and angularity unmeasured in the fourth (k = i or j) and fifth 
(k # i, /) rows. In the second-to-last row we summed the first two rows over all 
pairs of jets i, j to obtain the measured contribution for a specific t*, and in the 
last row, we summed all unmeasured soft gluon contributions. In the last two 
rows, we have taken the large t limit, jp = 1 in all cases. 



The non-cusp parts of the anomalous dimension of 5' meas and 
£ unmeas gjj^g t jj e same dependence on r, and therefore we are 
free to shift non-cusp terms freely between anomalous dimen- 
sions. While this does not change the physics, it allows us to 
organize the anomalous dimensions to match the contributions 
in Table [T] which we find more convenient for assembling the 
solution to the soft RGE Eq. (29) . By making the non-cusp 
part of 5 meas (r*) zero, we find that the shifted S meas (r*) is equal 
to SHii) from Table [l] and that the shifted S unmeas is equal to 

Finally, we can give the soft function anomalous dimension. 
Omitting terms which are suppressed by 0(1 It 1 ), the soft func- 
tion anomalous dimension is 



7s(r l a , 
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The solution of the RGE is 



S(t u ...,t m ;h) = Jdr\ ■■■dT' M S(T\,...,T , M ;jj Q ) 

M N 

x t/ s pai V, W) )f[ U k s (n -r' k ;fi,no) f[ U^,^), 



(32) 



!=M+1 



where Ug(jn) is an evolution kernel of a convoluted RGE and is 



and U' s and U^ dlr are evolution kernels 



of the form in Eq. ( fl9) 
of multiplicative RGEs and are of the form in Eq. ( |15) , The 
evolution kernels U^(Tk), U' s , and t/£ air correspond to the soft 
anomalous dimensions from S*(t*), S', and S pair in Table [I] 



5.4. Consistency of Factorization 

Adding together all jet and soft anomalous dimensions, we 
find, miraculously, the R dependence cancels between the un- 
measured jet anomalous dimension Eq. ( |25] > and sum over un- 
measured jets in the soft function Eq. pT| ), and the t„ + de- 
pendence cancels between the measured jet anomalous dimen- 
sion Eq. ( p6) and the sum over measured jets in the soft func- 
tion. The remaining pieces precisely match the hard anoma- 
lous dimension y# given in Sec. |5.1| such that the consistency 
condition Eq. ([3]l is satisfied. Note, however, that satisfying 
Eq. Q exactly required that we drop corrections of 0(1 ft 2 ) 
in the soft function. Requiring consistency of the anomalous 
dimensions at one loop has provided the measure t 2 » 1 to 
quantify the condition we used in justifying the factorization 
theorem in Sec. [3]that jets be "well separated". 

6. Application: Jet Shapes in e + e~ -» 3 Jets 

As an example of using the above results to calculate a jet 
observable in an exclusive multijet final state, we give the re- 
summed angularity jet shape distribution for a single measured 
quark or gluon jet in a three-jet final state in e + e~ annihila- 
tion. The techniques to derive and solve the RGEs to resum 
logarithms in jet shape distributions in SCET are essentially 
identical to those for event shape distributions as performed in 
US [35] EE 

We assemble the appropriate RG-evolved hard function, mea- 
sured jet function, two unmeasured jet functions, and soft func- 
tion given in Secs.|4]and|5] Evolving these from their tree-level 
values at initial scales ^h,^'j,^s to the scale /j. with NLL run- 
ning, we obtain the distribution in the shape r a of jet 1 with jets 
2,3 unmeasured, 



1 do- Pl p 2 p } 
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where crp^Pj is the cross section differential in the three jet 

-p2 j2 J_ 

momentaP; = a>;n,-, the effective hardscalewn = (o) x 1 w 2 2 oj 3 3 )t 2 
where T 2 = T 2 + T 2 + T 2 , and 'K is the sum of the hard, jet, and 
soft evolution factors, 



% = KhQx^h) + YjiK^fy + K^fis)] + Kf r (^Vs) ■ 

(34) 



1=1 



Inspection of Eq. ([33) suggests the reasonable choices for initial 
scales 

(i H = 0) H , fij = uiT l J a ~ a \ fx 2 / = w 2; 3 tan -, fx s = (DiT a . 

(35) 

For the unmeasured jet scales /j. 2 j 3 we kept in mind the factor 
of In tan 2 | present in K 2 , (see Table [TJ. To obtain the shape of 
a quark or gluon jet from Eq. ([33) we designate jet 1 as either 
quark or gluon and plug in the appropriate color factors and 
anomalous dimensions from Table[T]into u>f and Kf appearing 
in Eq. ( [33) . We report on a more detailed phenomenological 
study of these jet shapes in 11251 and their application to the 
discrimination of quark vs. gluon jets in future work. 

7. Conclusions 

We have demonstrated the intricate fashion in which the 
factorized cross section to produce exclusive N-jet final states 
when M < N are measured with a jet observable remains con- 
sistent for NLL running. We identified sources of power cor- 
rections to this factorization theorem and the consistency con- 
dition. Up to these corrections, the factorization theorem re- 
mains consistent independently of the number of measured and 
unmeasured jets and number of quark and gluon jets. 

One novel power correction that explicitly manifested itself 
in our calculation is in the separation parameter t . Since 1 jt is 
identically zero for all jet sizes when jets are back-to-back, this 
parameter has not been identified in the literature before. 

We find that, when a jet measurement is performed, the 
NLL resummed result has no dependence on the jet algorithm 
across the algorithms we considered (the Snowmass and SIS- 
Cone cone algorithms and the inclusive kx, anti-kx, and the 
Cambridge-Aachen kx-type algorithms). In addition, for un- 
measured jets the dependence on the jet algorithm parameter R 
(or D) is universal across these algorithms at NLL. 

Jet shapes such as angularities can be used to describe the 
substructure of a jet, and can be used, for instance, to distin- 
guish quark jets from gluon jets. In a future publication we will 
develop and describe a strategy to do so. We presented our cal- 
culations in such a way that allows for straightforward adapta- 
tion to other measurements as well, as we separated those parts 
of the jet and soft function that depend only on the jet algorithm 
and not the choice of jet observable. In addition, the ideas we 
discussed such as the power corrections that arise in the fac- 
torization formula and the method of calculating the soft and 
jet functions, will carry over to a calculation involving jet al- 
gorithms at hadron colliders, essentially amounting to having 
algorithm parameters that are invariant under boosts along the 
beam axis. 
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